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AH INVESTIGATION 07 AIRCRAFT HEATERS 
XIV - AN AIR AND HEAT PLOW ANALYSIS OP A 
RAM— OPERATED HEATER AND DUCT SYSTEM 
By L. U. K. Boeltor, E. H. Morrin, R. C. MartinelH, 

and H. P. Poppendiek 



SUMMARY 



A method cf graphical analysis is outlined which 
allows the prediction of the thermal and aerodynamic per- 
formance of a ram— operated heater and duct system, for 
cabin heating or wing de-icing, at any altitude and air- 
plane speed. This performance may he predicted from the 
isothermal tptal pressure drop across the components of 
the duct system at t uttsI 1 rent ile/t ing air r*t»«,:anfl th* 
thermal output of the heater at various air and exhaust 
gas rates. 



INTRODUCTION 



The performance of a. ram— operated aircraft heater is 
determined largely hy the amount of ventilating air which 
can he forced' through the heater— duct system by the ram 
pressure* The ventilating air rate, at a fixed airplane 
speed and altitude, depends upon the resistance to the 
flow of ventilating air, not only of the heater hut of 
the complete duct system from air scoop to the point of 
final air discharge. ■ Tor a given duct system, the pres- 
sure drop at a fixed air rate depends on the density cf 
the air passing through the system, and thus, at any given 
altitude the pressure drop for a fixed ventilating air 
rate is determined by the temperature of the air leaving 
the hoater. Conversely, for a given ram pressure, the 
resultant air rate through the duct system depends upon 
the temperature of the air leaving the heater. Because 
the temperature of the air leaving the heat exchanger at 
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a particular vontilating air rate dopen&B upon tho thermal 
output of the heater, the determination of tho per f or mane o 
of a heator duct system involves the simultaneous solution 
of a prossuro drop equation and a heat transfer equation, 

A graphical solution of these equations, for tho de- 
termination of the weight rate of air passing through the 
heater—duct system and the temperature of the air leaving 
the heater, at any altitude and airplane speed, is pre- 
sented in this report. If the desired ventilating air 
weight rate at a certain altitude and plane speed 1b known, 
the allowable duct losses may be determined by the method 
presented. 

It is well recognized that, in many cases, the limi- 
tation on the design of a her.t 'exchanger duct system is 
the question of allowable space and so forth, and in Buch 
cases the analysis presented below may not be pertinent, 
although it may be used as' a guide for design. 

Shis Investigation, conducted at the University of 
California, was sponsored by, and conducted with financial 
assistance from, the National Advisory Committee for Aero- 
nautics . 



SYMBOLS 



A cress-sectional area of flow, ft 3 

A a cross— sectional area of flow at section a, ft 3 

At, cross— s ect ional area of flow at section b, ft 3 

Ab. constant cr os s— sect ional area of flow in heat 

exchanger, ft a 

A 8 oross— sectional area of flow at section 2—2 at 

inlet to air scoop, ft a 

A 3 cross— sect ional area of flow at section 3—3, 

entrance to heat exchanger, ft 3 

A 4 cr osb— sect ional area of flow at section 4—4, 

exit from heat exchanger, ft 3 
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A B cross— sect lone, 1 area at section 5—5 at entrance to 

isothermal discharge section, ft 3 

£ 8 cross— sectional area of duct at Bectlon 6—6 at 

point of final air discharge, ft« 

c p heat capacit7 of air, Btu/lb °F 

D hydraulic diameter of duct, ft 

F frictional pressure loss, lb/ft a 

k Btu 1000 3tu (kilo Btu) 

g gravitational force per unit mass, 32.3 Id/ (lb sec a /ft) 

L distance along duct, ft 

n exponent obtained from £F ft b against V± BO curves 
(See references 3 and 5.) 

F absolute static pressure, lb/ft 2 

P a absolute static pressure in free air stream at any 

altitude, lb/ft* 

P_ absolute static pressure at inlet to air scoop, 

lb/ft s 

F 3 absolute static pressure at entrance to heat 

exchanger , lb/ft 5 * 

F 4 absolute static pressure at exit of heat exchanger, 

lb/ft a 

F B absolute static pressure at entrance to isothermal 

discharge section, lb/ft a 

P absolute static pressure at point of final air 

discharge, lb/ft 2 

F. mean absolute static pressure in duct system during 

80 the isothermal total pressure test (or that based 

on calculations using appropriate data), lb/ft 3 

p static pressure, lb/ft 8 

3-alt thermal output of heater under operating conditions, 
Btu/hr 



q lal} thermal output of heater as measured in laboratory, 
a Btu/hr 

B gas constant in PV = BT, ft lb/lb °F 

T absolute temperature, °E 

T a absolute temperature at section a, °E 

T-jj absolute temperature at section b, °B 

T £ 1 absolute tenperature of air in free air stream, °E 

T absolute temperature of air, just inside air scoop, °E 

t „ temperature of air at entrance to hsat exchanger, °F 

3 

T 3 absoluta temperature of air at entrance to heat 
exchanger, °B 

t 4 temperature of air leaving heat exchanger, °P 

T 4 absolute temperature of air leaving heat exchanger, °B 

I g absolute temperature of air after p.'-.ssing through 
non-isothermal duct, at entrance to isothermal 
discharge section, "£ 

absolute temperature of air pi.spir.,3 through duct 
130 system during isothermal total-prr-ssuro test, n H 

u mean v-olocity of flow at any section of duct, ft/sec 
m 

u t velocity of air stream relative to airplane, ahead of 
air scoop (true airspeed of ainlane + air 
velocity produced by propeller), ft/ssc 

u a velocity of air relative to airplane at section 2-?, 
entrance to air scoop, ft/s.^c 

u 6 velocity of air relative to airplane at point of 
final air discharge, ft/sec 

V specific volume cf air (l/density), cu ft/lb 

V s-jecific volume of air at section a. cu ft/lb 
a 

specific volume of air at section b, cu ft/lb 

V -specific volume of air in free air stream, cu ft/lb 

V specific volume of air Just inside nir scoop, 

cu ft/ lb 



I 
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Y" 3 specific volume of air at entrance to heat 
exchanger, en ft/l"b 

V- specific volume of air at discharge from heater, 
cn ft /Id 

"V B specific volume of air after passing through non- 
isothermal duct, ou ft/lb 

V 8 specific volume of air at point of final air dis- 
charge, ou ft/lb 

T lao specific volume of air during isothermal total- 
pressure test, cu ft/lb 

W & ventilating air rate, l"b/hr 

tf exhaust £as rate, lb/hr 

V Jgo ventilating air rate during isothermal total- 
presTure teat, lb/hr 

x distance along duct, ft 

AL distance alonp duct, across which AP is 
measured, ft 

AP . isothermal frictional Tirrssura los3 betwoon section 
a a and section b, lb/ft a 

AP a frictional pressure loss "between freo air stream 

p.nd entrance of n,ir scoot?, for isothermal condi- 
tions specified V/ P ieo ; T lao , V t _ eQ , lb/ft a 

AP g frictional pressure loss "between entrance of air 

scoop &nd entrance to heat exchanger, for 

isothermal conditions specified "by P-i,,... E t , 
W iE0 , lb/ft a 

AP 3 ^ 4 frictional pressure loss across heat exchanger for 

isothermal conditions specified "by Pi an » T iio» 
and W iao , lb/ft 8 

AP 4 _ g frictional pressnro loss along the air duct fro^i 

the heat exchanger to the entrance of isothermal 
dischargo section for the isothermal conditions 
specified by P ia0 , T iao , and tf iao , lb/ft a . 
The isothermal prassure loss AP 4 __ B is made 

up of the sum- of the pressure losses through all 
bends, sudden expansions, straight duct, rnd 
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so forth, which occur along the duct from the 
heater to the entrance to the is othermal. dis- 
charge section. - 

&F S _ S frictional pressure loss in isothermal discharge 
section for the isothermal condition specified 
*y p iso> T iso> and * le0 . lb/ft 8 

AL V-m ^ 

t isothermal friction factor defined by A3? = £ — — — — 

D 2g 

T a ^ temperature of air entering heat exchanger, °3? 

T„ temperature of hot exhaust gases entering heat 

exchanger , 3P 

K. Q frictional pressure loss coefficient for fluid 

expansion or contraction defined "by equation (5) 

K- 0 frictional pressure loss coefficient for bends , 

and so forth, defined by equation (5) 

k exponent for ad,iabatic compression defined by 

equation FV k i P i V g k 



BASIC EQTJAIIOiTS 



5?he derivation of the equation for the prediction of 
the non— is other mal pressure drop through a heater— duct 
system, follows directly from the integration of the 
modified Bernoulli equation along the flow system. 



The basic elements of a heater— duct system are shown 
in figure 1, total— head tubes being visualized as located 
at sections 1 to 6. 

two distinct types of flow system are indicated in 
figure 1: 

Slype I *-<• Isothermal flow — unequal cr oss- sectional 
areas at points of total pressure 
measurement 



2ype II 



Hon— is other mal flow — unequal cr oss— secti onal 
areas at points of total pressure measure- 
ment 
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The components of the flow system shown in figure 1 
may he classified as follows: 

1. The approach section, in which the air is adia— 
batically compressed hy the motion of the airplane through 
the air. The flow is non-isothermal. 

2. The inlet section, in which the cross— sectional 
areas at the two total head stations 1 and 2 are different 
and the flow is isothermal (type I) 

3. The heater, which in the general case is assumed 
to hare unequal cr os s— sect ional areas at the two total 
head sections 2 and 3, and in which the flow is non- 
isothermal (type Jl) 

4. The discharge duct, in which the cr oss— sectional 
areas at points 3 and 4 differ and the flow is non— iso- 
thermal because of heat losses along the duct (type II) 

5. The final discharge seotion is of the same type as 
the inlet section (type I) 

The Bernoulli equation in differential form (including 
friction, hut neglecting elevation differences) is: 

- V dp - -— + V 41 (1) 
6 

A second convenient method of writing equation (l) is: 

. d p = (-JL-) JL ^ d 'I s - + djr (2) 
p KasOoJ gT W \As 

Equation (2) results from the substitution of tL. » — ~r-r 

^ m 3600 A 

into equation (l). The significance of the terms in equations 
(l) and (2) is as follows: 

1. The first term represents the difforenoe in static 
'prassvi.ro due to both the changes in kinetic energy of tho 
fluid and frictional losses. 

2. Tho second term represents tho change in prossuro 
duo to changes in the kinotie energy of the fluid, which 
may result from: 
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(a) Changes in density or specific volume (V) 

(b) Changes in cross— sectional area of flow (A) 

3. The last term represents the change in pressure 
due to frictional losses. The latter Iobb is irrecoverable , 
since the energy lost by friction is dissipated in the form 
of heat, 

liquations (l) and (2) may be integrated between any 
points in the flow system subject to the restrictions 
imposed upon the use of the Bernoulli equation. The in- 
tegration of the second term in equation (2), in particu- 
lar, yields different results for the two types of flow 
system discussed. 

T yp e I. Iso t herm a l flow — un equal cross— sectiona l 
ar eas a t points of t o tal -pr e ssure mea s urement . — 

f ( JL.Y i fIS 4 CSS . :'-2L_Y J- - (2>T] W 

J gV \kj W \3600> 2g LU b / \A J . 

a 

Type II. Son— is othermal flo w — unequal cross— sectional 

a reas at point s of t otal pressur e m eas u rement .— J or this 

case the integral cannot be evaluated precisely, since the 
integrand is not an exact differential as it was for 
case I. Several methods of approximation may be used. 

(a) Average spec i fic v olume .— If the flow system is 
near ly isothermal, the following approximation may be used: 

i } ( v Y i N\ . (i\ ( w V i iVv^ 8 ^yi 

a 

(b) Average c ro ss— sect ional area .— If the area is 
almo st constant, the following approximation may be used: 

_i_ } ( v n b 1 4 .'A i 7I£ _ fV\l 

A av / ^3600/ g \LJ ~ V,3600; g A av \±\J 

(c) If the area is constant, no approximation is 
necessary, for the integral can be evaluated exactly. 
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a 

(d) Tor the system shown "below, which represents a 
heat exchanger of constant cross- sectional area A^, 



A a 


r 


Ah 


Ah 


i 






T a 









application of equations (5) and (4c) yields: 



r ( t v i ^ y v 



[(ft* 



An isothermal contraction (or expansion) at point a, 
a non— is o thermal change of the fluid at constant cross- 
sectional area, and then an isothermal expansion (or con- 
traction) at point b is considered for the derivation of 
the preceding expression. 

The third term in equations (l) and (2) may he 
evaluated as follows: The irrecoverable frictional loss 
for flow in a conduit usually consists of: 

1. Skin friction losses in straight sections of the 
duct 



2. Losses due to sudden expansions and contractions 
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3. Losses due to flow around "bends , and so forth 
The skin friction loss usually is expressed by 

d3 p = £ ?s! £Lx = I Y 1- dx 

1 2g DV \36OO / 2g DA 

where £, the friction factor, is a function of the 
Reynolds nv.mber of the flow system. 

Th 1 lossos sudden expansion pnd contraction 
usually are expressed by an eq.uetion of the form 



3 

dF a = E _r?_ = K 

9 2£ T 6 V 3600 J 2- A 3 



where the expansion coefficient Z e m».y bo estimated 
from data in references 1 and 2. 

The losses fue to flow around bends r.lso are ?xpress:sd 
by an equation of the form 



dr- = K-u. --S- = Z 
3 T 2s 



2a! - z f __L V JL, 
is 13 \ 3600 / 2* A* 



Thus, for a duct which includes several banc's, sudden 
expansions, and 30 forth (see references 1 *nd 2), 



dF 



= (-7-- ^ 3 I ^ --o dx + ... + K -I + ... J- + ...-I. 

V3600 / 2;t L DA 0 A a D A 3 J 

Int -3gratin,*» along the whole duct length for conditions of 
W isc' P i 3 o» T iso yields 

a ^ \36oo/ 2 C L--' DA 3 <^'A 3 <^ a^J 3 

A plot of A^ao against W i so usually will reveal an 
equation of the form 

av - v Via0 ( W iso \ n 
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where Z varies slightly with temperature "because of the 
variation of Reynolds number with temperature and the 
resultant change of the friction factor £. The exponent 
n will vary somewhere between 1.75 and 2.00 due to the 
variation of the friction factor i with V. If the 
greatest friction loss is due to expansions and contrac- 
tions, n will be near 2. If the greatest loss is due to 
skin friction, the exponent will be near 1.75. Typical 
values of n are given in references 3 and 4. 

In reference 3 it has been shown that the isothermal 
pressure loss AT a b can 06 corrected to non— isothermal 
conditions by correcting for changes in specific volume 
and viscosity with temperature. If the exponent n is 
known, the non— is othermal frictional pressure loss at any 
weight rate W is therefore: 

/ v \ a /T +T \°« 13 /y +v \ 



isothermal frictional pressure loss between section 
a and section b at the temperature ™i B0 i air 

specific volume v i 8 o» and weight rate w i 0O » 
lb/ft* 

weight rate of air, lb/hr 

exponent (between 1.75 and 2.0) obtained from plot 
of AP al) against tf la0 

arithmetic average temperature in length of duct ab 

arithmetic average specific volume of air in length 
of duct ab , cu ft/ lb 



where 
APab 



n 

T +T. 
a ] 

2 
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INTEGRAT IOiJ OT BERNOULLI EQUATION ALOirO HEAT EE- DUCT STSTEIT 
Pros Stream to Entrance of Scoop 
(Sections 1 and 2) 

PI ow_ s jfa t em_-_(Juneo fr - : ial_floT.'_ areas i_co.5?;2res sible _f^2?l • 

In order to show the effect of the compressibility of the 
fluid betveen sections 1 and 2 the frictionp.l loss AF 

l— e 

is postulated t - be ne^l i£ibl e . Compressibilit - ' - effects, 
ovi the other hand, c-mnot he neglected. Equation (l) then 

h 3C07.83 : 

u_ du_ 
- V dp = -E___S 

e 

' The thermodynamic path followed "by the i»ir passing 
from section (l) to section (2) is postulated t -> he 
adiabatic. Then: 



= P^* 



Thus 



P 3 ;j p 8 _ £-i 

r -Y d-n = - P k T. / S? = p T 1 1 - fell" £ 




and from equation (l) 

k-i 



k 



■:~ v«L l " <^;) J = -5c" 1 - 



(7) 



An approximate solution of the preceding equation for 
(? 3 - ^3^)1 which is accv.rate to within 10 percent if 

P - P 

— — < o.H is presented below: 
P 

1 

Because 

la = J?a-:.!i2 = 1 + 

P i P ! P a 
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Koglootiag all but the first two torms, and subetitut ing 
- in equation (7) yields 

- (P a - Pi ) , ^ i =~^ (8a) 

which is the solution that would result by considering the 
fluid incompressible. 

The frlctional pressure loss A3T 1-a (which was neg- 
lected in determining the effect of fluid compressibility) 
is added to equation (3a) for use in the following flow 
.equations, yielding 

Although the pressure P a is not greatly influenced by 
compressibility effects, the temperature at point (2) 
may be several degrees higher than that of the free air 
stream due to compressibility. She temperature T a ma;' 
be readily estimated. 

Substitution of the relations 





= RT 1 










• v k 
1 1 





k 

s 



in equation (?) yields 

u a a - -a-i 3 - H k T 1 
3g k - 1 ta 1 a J 



then 




or 
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*. - «» ♦ W c.) 

Scoop to Heater 

(Sections 2 and 3) 

Plow system — type I ( I sotherma l , uneq u al flow areas ) • — 
Integrating equation (l) "between sections 2 and 3 gives 

or, since u = g60Q A and T fi = T 3 , V E = V 3 

% -* -^OTJ 2iXI " 2gTI + ^ a "VW \jiTo) VW (11) 

Adding equations (ll) and (8b) in-order to eliminate u a 
results in the equation* 

V^iao/ N^loc^ \Jisoy 
u, a 

ITote that P + 9 * . is equal to total pressure in 

1 og» j 

free air stream "before scoop. 

Heater 

(Sections 3 and 4) 

Tiow syst e m — type II (ri on— is othe rm al . un equal flow 
a reas ) .— The integration of equation (3) between points 3 
and 4 fields : 



*A term involving the difference in the specific volumes 
at sections 1 and 2 is neglected in the determination of 
equation (12). 
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*3 



-'•-(sfc)W®*0S-G5-')J 

•■M&jseras) 



Discharge Duct 

(Sections U and 5) 

Flow system - type II (non-iaothernal uneq.ual_f low 
arfiflfl).- Integration of equation (2) yields 

Final Discharge Section 

(Sections 5 and 6) 

Tlo'rf system -_type I ( is otherraal , unequalar eas ) 
Integration of equation yields 

. _ 8 « .» /If \ n / T«t V' 13 / Y« \ 

P. 



5 



*• " r g - T ; " sr. + ATb - Cr~) fe) fey (15) 



or 

a 



?B " ( Pe + asYg)" '".(3600/' igI7 . 

where P 6 + — is equal to total pressure at point of 
air discharge. 
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Eliminating the intermediate pressures P 3 , P 4 , and 
P B between equations (12), (13), (14) , and (15) results 
in the equation for the pressure drop "between points 1 
and 6 

( P i + 5e0 " ( Ps + ^0 = (wot) A 

[©■ - cai - a * v-ae;-"Ge) 

By. the use of the ideal gas equation PV = ET equation 
(17) may "be simplified 




(18) 
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In equation (18) the terms on the left of the equal 
sign represent the difference in total pressure "between 
the free air stream and the point of air discharge. The 
first term on the right of the equal sign represents the 
pressure changes due to the acceleration of the air in . 
the duct, which are due "both to changes in area and changes 
in specific volume The last term represents the irrecov- 
erable pressure loss due to the friction in the complete 
duct systems. It should be noted that each isothermal 
frictional pressure Iobs is corrected to the' operating 
temperature by different temperature corrections, 
depending on the type of flow system represented by each 
separate AT. Thus, any complex flow system cab he 
broken up into a series of systems, and the pressure 
drop through eaoh corrected to n on—isothermal condition 
by the method outlined. 

In equation (18), for a given duct system for which 
the isothermal total pressure dropB AT la . a , AT a _ 3 , AT 3 __ 4 

AT 4 _ B , and AT 5 _ a are known, the remaining unknowns are 
W and T 4 • The fixing of the altitude, the airplane 
speed, and the heat loss from the duct establishes all 
other variables in the equation. Thus, for any altitude 
and airplane speed a curve of V against T 4 can be 
drawn which will reveal the rate of flow possible through 
the duct system for any temperature T 4 . 

2he relative importance of the various portions of 
the duct system may be readily established, for the largest 
of the corrected pressure drop terms in equation (l8)- r will 
be the term which controls the rate of air flow. If it 
becomes necessary to increase the rate of flow through the 
heater— duot system, attention should be focused on the 
largest term. By breaking up a complex duct system into 
a series of small units, the units causing difficulty 
then may be readily isolated. 

« 

SoJtas In many cases, the first sum of terms on the 
right-hand side of equation (18) (which involves pressure 
drops due to changes in area and fluid temperature) may 
be neglected when compared with the second sum of terms 
involving the frictional pressure losses (AT). 

After the curve of W against T 4 is established 

from a consideration of the pressure drop characteristics 
of the duct system (from equation (18)), the thermal per- 
formance of the heater must be utilised in order to 
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establish the operating point of the heater— duct system. 
This thermal performance of the heater then is used, to 
establish a second curve of V against T 4 , which is 
fixed by the thermal output of the heater , since for any 
particular W and ' eihaus t— gas temperature only one 
ma&nitude of T 4 is possible. The relation 



or 

= + * 3 (2°) 

WCp 

is utilized to obtain this second curvp. The heater 
capacity qi a b • determined in the laboratory, must 

be corrected to altitude and temperature conditions 
by the method outlined in reference 5. The intersection 
of the ourve of W against T 4 obtained from' the pres- 
sure drop characteristics of the heater— duct system 
(equation (18)) and the cur-re of W against T 4 from 
equation (20) fixes the operating point of the system 
at the particular altitude and airplane speed under 
consideration. A series of calculations at various 
altitudes and airplane speeds then will establish the 
complete performance of the unit. 



Eecapitulat ion 

The performance of a ram— operated heater— duct system 
depends on the pressure loss characteristics of the complete 
system p.nd the thernal output of the heater. 

The following data must be known in order to estab- 
lish the performance of the system: 

1. Thermal output of the heater as a function of air 
rate at various exhaust gas rates and exhaust 
gas temperatures. The methods of correction 
of the performance of the heater to any alt i^^:'.-- 
tude are presented in reference 5. 
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3. Isothermal total pressure loss through the vari- 

ous pertinent portions of the air duct system. 
Total pressure data are necessary in order to 
evaluate the irrecoverable frictional pressure 
loss . 

3* Airplane speed, altitude, air temperature, and 
static pressure in the free air stream 

4. Heat loss from discharge duct 1 " 

5. stftt^c preBsure at the point of flnai air discharge 

If these quantities are fcnown , the complete perform- 
ance of the heater— duct system may be calculated at any 
altitude. Or, if the weight rate V is fixed by design 
at a .certain altitude and airplane speed, the allowable 
isothermal total pressure loss for any section of the 
duct system may be calculated. 



EIA.HPLB 



An example of the application of equation (18) to 
the prediction of the performance of a ram— operated 
exhaust gas to air cabin heater and associated duct woric 
as t. function of airplane speed and altitude If presented 
balov: 

Data 

The following data are available: 

1 . Therm al output of heater (fi g. 2).— The data 
shown In figure 3 were obtained in the laboratory for the 
following conditions: 

(a) Atmospheric pressure, 14.7 ps ia 

^ / \ o 

! (b) Temperature of air entering heater, T a ; 100 T 

1 " 1 

(c) Temperatures of exhaust gas entering heater, 
T c , 1400° JP 

U f i_ _ ; , 

' *The duct loss oan be estimated by means of the equations 

presented in reference 3. 



L 
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Tor purposes of the analysis the following flight 
conditions will "be utilised:* 



Altitude 

10,000 
20,000 
30,000 
40,000 



1600 
1600 
1600 
iGoo 



+ 23 
-12 

-48 
-67 



\ 

4000 

Uooo 

4000 
4000 



With this data and the data in figure 2 the heater output, 
and thus the temperature of the air leaving the heater, 



4* 



may he calculated readily for any air rate, » a , 



through the heater. This procedure will yisid one of the 
curves of W ft against t 4 required in the analysis. 

2. The isothermal total pressure drop data for the 
heater-duct system for an air rate of 30 n 0 -Hounds per 
hour are shown in figure 3» ^ ne exponent n = 1.8. 

3. The performance of the unit is desired it the 
following airplane speeds and altitudes: 



Alt itude 



(ft) (deg H) (lb/ft 3 ) (cu ft/lb) 



Airplane speed 
(true air speed) 

LW2*l 



10,000 
20,000 
30,000 
40,000 



483 
448 
412 

393 



14S0 

972 
627 

392 



17.7 

24.5 
35.0 
53.5 



100 
200 
300 
400 



4. Heat loss from discharge duct is 1000** Btu/hr ft a . 
Surface area of duct is 2^ ft 2 . 

*Tor lack of better information, and in order to reduce 
the complexity of the sample calculation, constant magni- 
tudes of Tg and 'tf g were used for all airplane speeds 



and altitudes. If data for the variation of t 



and 



Wg with altitude and airplane speed are available, they 
may be readily utilized in the analysis in place of the 
constant values employed in this report. 

**The duct loss can be estimatod by means of the equations 
presented in reference 3* 
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5« Static pressure in cabin is equal to atmospheric 
pressure at the given altitude. 



Sample Calculation 

A sample calculation is performed for an altitude 
of 30,000 feet. 

1. Calculations of weight rate V as a function 
of t 4 from pressure drop data at an airplane speed of 
3Q0 miles per hour. Equation (IS) is utilized. The 
various terms in equation (IS) have "been assigned the 
following magnitudes for this calculation: 



Sour c e. 







P e » 627 lb/sq 


ft (U.35 lo/sq in. ) 


Design condition 




= 


300 mph = UUl ft/sec 


n it 


T l 


a 


35 cu ft /lb 




it n 


■g 




32.2 ft/sec a 




it n 


H 


a 


53.3 ft 




n n 


n 


e 


1.8 




n n 


*i 




Ul2° R 




n n 




a 


T 3 = U27 0 E 




Equation (9) 


A3 


B 

=1 


0.235 ft 3 
0.225 ft 8 




Figure 3 

Tf It 


A 4 




0.2^5 ft 3 




n n 


A B 




O.I96 ft a 




n 11 


±e 




0.500 ft a 




it 11 




IS 


_ _ 2U X 

8 4 0.2U 


1000 
x V 


Design condition 


C P 




0.2^ Btu/lb °3P 




Reference 5 
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Somas 



W iso 


= 


3000 


lb/hr 




Figure 


3 


*i 8 0 


= 


532° 


B, 




it 


n 


P iao 




2120 


lb/ft 3 


= 1U7 lb/sq in. 


it 

■1 


11 




B 


0.52 


lb/ft s 


(0.1 in. H a 0) 


n 


it 






0.5?. 


lb/ft a 


(0.1 in. E a 0) 


n 


11 






6.23 


lb/ft a 


(1.2 in. E a 0) 


n 


it 






11.7 


lb/f t a 


(2.25 in. E a 0) 


11 


11 




3 


2.69 


lb/ft a 


(0.50 in. E a 0) 


n 


it 



Substitution of theae quantities in equation (lg) 
yields an equation in V (the ver.ti lat in~ air rate) 
and T 4 (the temperature of the air leaving the heater). 
Chooaing arbitrary values of t 4 = 200°, U00°, 6C0°, 800° 
and 1000° I" results in the following aquation in W. 
The coefficients A and B are the multipliers of W a 
and V 1 * 8 (since n a l.g), r es-necti vel~, in equation (lg) 

S6.5 «= A \I a + B V 1 ' 8 

where 



t 4 A x 10 8 3 x 10 8 V 

L^BS-ll (lt/hr) 

Altitude 200 6l U510 3070 

30,000 ft Uoo 105 5920 2620 

and 6C0 1U9 7360 232O 

Airplane speed, 300 205 3790 2110 

300 mph 10C0 237 10300 I930 



A plot of the curve against t 4 obtained in the 

foregoing is shown in figure together with the curves 
calculated for airplane speeds of 100, 200, and h00 miles 
per hour. 

The heater output must now be utilized to obtain the 
operating points of the h-.at er-duct system. From figure 
2 the following data are obtained: 
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Wo 

(lb/hr ) 


^lab 
(Btu/hr) 


^alt 
(Btu/hr) 


*alt 
Ve P 


t . 

4 

(de fi J) 


1000 


150,000 


188,000 


783 


751 


2000 


215 ,000 


270,000 


562 


529 


3000 


27 0,000 


340,000 


472 


43 9 


4000 


310,000 


390,000 


406 


372 


5000 


337,000 


424,000 


363 


.317 



Altitude , 

30,000 ft. 

and 

Airplane speed, 
300 mph 

The laboratory data were corrected to altitude by multi- 
plying by the ratio:* 

/ 1600 -t s \ 
4 alt = q lab U400 - XOoJ 

The temperature t 4 was calculated from the relation: 



l alt 



W c-, 



<n - 



(21) 



The resulting curve of W against t 4 is plotted in 

figure 4. The intersection of the curve obtained from 
the heat transfer data** and the curves obtainod from 
the pressure drop data yields the operating points for 
the heater— duct system at 30,000 feet altitude. This " 
performance is summarized as follows: 



•This correction neglects the small changes in heater out- 
put due to changes in the unit conductances on the air and 
gas sides resulting from the fact that the average oper- 
ating temperatures of the air and exhaust gases are dif- 
ferent from those utilized in the laboratory test (referen 
5). The tomperaljure t 3 is equal to the outside air 
temperature (—48 T) plus the temperature rise due to 
adlabatio compressibility. 

**If variations of heater output with airplane speed due 
to changes in exhaust gas rate and temperature are to be 
includod in the computations, a separate curve of V 
against t 3 from the heator output data will be obtainod 

at each airplane speed. Since the exhaust gas rate, is 
assumed to be constant in these calculations, tho hoator 
performance is represented by a single curve in f ig f 4. 
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Altitude, 30,000 ft 

Alrplr.no speed W t 4 ^alt 
(m ph) ( lb/hr) (d og J) (Btu/hr) 

100 550 960 133,000 

200 ■ 1420 625 230,000 

300 2450 475 310,000 

400 3600 3S0 380,000 

Tho repetition of the preceding procedure for alti- 
tudes of 10,000, 20,000, and 40,000 foot yiolds the com- 
plote porformance of tho unit. The final results are 
tabulated in table I and plotted in figures 5, 6, and 7, 



DISCUSSION OF RESULTS 



figure 5 reveals the temperature at which the "air 
will lo.'.vo the heater at various airplane speeds and 
altitudes* Inspection of tho figure shows that the lower 
the airplane speed the higher the temperature of the air 
loaving tho heater, and tho greator tho altitude, the 
higher tho tomporature leaving the heater. All the cal- 
culations for figure 5 wore basod on a vontilating air 
duct in which no air— control vpIvob woro installed, and 
thus figure 5 represents tho optimum pprf ornance of tho 
unit. Installation of valves for restricting tho rate of 
air flow will, of course, have tho same effect as in- 
crossing the rosistanco to flow of tho duct.* 

figure 5 represents the air rate through the heater— 
duct systoui at various altitudes and air speeds. As 
would bo oxpected, the lower the altitude and tho higher the 
■riir^lane speed,, tho greater tho air rate. Combination 
of the data prosonted in figures 5 and 6 allows tho pre- 
diction of hoator output at various altitudes and airplano 
spoods. Thoso results aro shown in figuro 7. 

Tho data in figures 5, 6, and 7 establishes tho por- 
formance of the hoator— duct system. The snmo calculations 
may bo readily porf or mod for a wing do— icing systom, the 

♦If a blower is installod in the vontilating air duct on 
the upstroam Bide of tho heater to supplomont the ram pros— 
euro, a term APf anL , equal to the pressuro head provided 
by the blowor, may bo added to tho loft side of equation 
(18). The romaindor of tho analysis romains unchanged. 
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main difference in the calculation "being a greatly lowered 
T B due to the heat lost "by the air ae it flows in the wing 
de-icing ducts ("between sections and 5). 

Oases "may arise in which the duct system is fixed and 
the performance of several heaters in this given duct 
system is to "be oompared. Repetition of the calculation 
shown for the various heaters will allow a rational com- 
parison of the performance of the several heaters to "be 
made* 

In other cases the heater output at a given alti- 
tude and airplane speed may "be fixed "by design, and the 
allowable duct losses must he determined. The following 
procedure may "be followed in this case: 

(a) From the known. heater output as a function of 
air rate, gas rates, and temperatures, the 
desired air rate V and air outlet tempera- 
ture t 4 may "be calculated "by employing 
equation (20). 

(h) Substitution of these values of V and t 4 

into equation (IS) allows the evaluation of 
one unknown isothermal pressure drop, say 
£iI' 4 _ B . The duct then must "be designed to 

function within this allowable -nressure drop. 



C0I:CLU3I0iTS 



A method has "been presented in this report for the 
prediction of the approximate thermal and aerodynamic 
performance of a ram-operated heater and duct system at 
any altitude and airplane speed, provided the following 
data are known: 



(a) Thermal output of heater at various air ratos and 

exhaust-sas rates 

(b) Isothermal total pressure drop (frictional Iosb) 

through duct system at several air rates 



(c) Heat loss from discharge duot (or win*, in 
de-icing system) 
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(d) Static pressure at the point of final air discharge 



University of California, 

Borkoley, Calif., December 1943. 
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TAB LI I 



Alt itude 

**■ JU V X U UUw 

(ft) 


Tr lift ftj.3TBT3ftftd 

of airplane 
(mph) 


t 

(deg ') 


v 

M J 

(lo/hr) 


q alt 
(Btu/hr) 


10,000 


100 


635 


1560 


228,000 




200 


430 


3750 


366,000 




300 


300 


6400 


426 ,000 




400 








20,000 


100 


785 


930 


179,000 




200 


500 


2400 


296 ,000 




300 


390 


4000 


386 ,000 




400 


300 


5850 


439,000 


30,000 


100 


960 


550 


133,000 




200 


625 


1420 


230,000 




300 


475 


2450 


308,000 




400 


390 


3600 


380,000 


40,000 


100 


1200 


320 


97,000 




200 


800 


800 


167,000 




300 


635 


1360 


230,000 




400 


510 


2000 


277,000 



Pi 



Air 
velocity 
u l 



P2 



P3 



P4 



P5 



- — — — €i 



Approach 

section 
Adiabatic 

compression 
Non-Isothermal 
Ai f Ag 



^Total headj tubes 




Note : 



Difference in isothermal total pressure 
equals frictional pressure loss, AF. 



/ 



P6 



~~1D 

Air velocity 



Note 



Inlet 

section 
Isothermal 
A 2 t A 3 



Heater 
Non -Isothermal 
A3 £ A 4 



Outlet duct 
Non-Isothermal 
A4 ^ A5 



Final discharge 
section 
Isothermal 
^ A6 



Figure 1.- Elements of heater and duct system. 
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Figure 2.- Thermal output of heater from laboratory 

data. (The heater output may be estimated 
by meajiB of the aquations presented in reference 3, 
if laboratory facilities are not available.) 




Ah 
A3 
A4 
A5 

A 6 



Areas 

0.235 ft 2 

0.225 ft2 

0.245 ft 2 

0.196 ft 2 
0.500 ft 2 



Isothermal flow data 

T iso = 532° R 

p iso = 14 - 7 p SIA 
V is0 =13.58 cu ft/lb 
W iso = 3000 lb/hr 

^1-2 = °- in H 2° 
AF 2 _ 3 =0.1" HgO 

AP3-4 = 1-20" HgO 

A]?4_5 = 2.25" HgO 

AF5-6 = 0.50" HgO 
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-Represents a 
total head tube. 



Figure 3.- Isothermal total-pressure drop data for heater and duct system. (The isothermal total-pressure 

drop may he obtained in the laboratory or may be estimated by means of the data presented in 
reference 4.) 
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t 3 = temperature of air leaving heater, °F 

Figure 4.- Curves of ff vs'tg at an altitude of 30,000 ft from pressure drop data (equation 18.) 

and from thermal output of heater (equation 20). The intersections of the curves 
represent the operating points of the heater and duct system at 30,000 ft altitude. 
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Figure 5.- Temperature of the air leaving the heater at 
various airplane speeds and altitudes. 
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Figure 6.- Ventilating air rate through heater 
and duct system at various airplane 
speeds and altitudes. 
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Figure 7.- Thermal output of ram- operated heater as a func- 
tion of airplane speed and altitude. 



